Testing the Cosmic Censorship Conjecture with point particles: the effect of radiation reaction and 

the self-force 
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A classical thought-experiment to destroy black holes was envisaged by Wald in 1974: it consists of throwing 
particles with large angular momentum into an extremal black hole, checking whether their capture can over- 
spin the black hole past the extremal limit and create a naked singularity. Wald showed that in the test-particle 
limit, particles that would be otherwise capable of producing naked singularities are simply scattered. Recently 
Jacobson and Sotiriou showed that if one considers instead a black hole that is almost, but not exactly extremal, 
then in the absence of backreaction effects particle capture could indeed over-spin the spacetime above the 
Kerr limit. Here we analyze back-reaction effects and show that for some of the trajectories giving rise to 
naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is 
important, and seems to have the right sign to prevent the formation of naked singularities. 

PACS numbers: 04.70.Bw, 04.20.Dw 



I. INTRODUCTION 

The appearance of singularities as solutions of a theory's 
field equations often signals the breakdown of that theory in 
some regime. In classical electrodynamics, for instance, the 
field of a point-like charge (the basic object of the theory) di- 
verges at the particle's location. This singularity can be re- 
solved only within quantum electrodynamics. Newton's the- 
ory also exhibits singularities at the location of point masses, 
and these singularities remain in General Relativity. In 1939, 
Oppenheimer and Snyder 1 1!] showed that under the assump- 
tion of perfect spherical symmetry, a sufficiently massive star 
that runs out of thermonuclear fuel will undergo indefinite col- 
lapse to a point, where the curvature invariants diverge. In 
fact, divergent curvature invariants as the outcome of collapse 
are a generic feature of General Relativity, even in the absence 
of spherical symmetry, as shown by Penrose in 1965 It 
is expected that a quantum theory of gravity will be able to 
smoothen and resolve this singular behavior 

Fortunately, the equations of General Relativity seem to 
contain a self-protection mechanism (from the unpredictabil- 
ity of naked singularities), hiding curvature singularities 
within event horizons and making them invisible to outside 
observers and therefore to experiments. This is known as the 
"Cosmic Censorship Conjecture", and was proposed by Pen- 
rose in 1969 m. While several counterexamples to this con- 
jecture are known, they either rely on using certain equations 
of state beyond their range of validity (such as the pressure- 
less matter equation of state used by Ref. iQl'); or they are 



non-generic in nature {i.e., they only happen for very spe- 
cific sets of initial data: see Ref. ijst]); or they are staged in 
higher-dimensional spacetimes Therefore, while a gen- 
eral proof of the Cosmic Censorship Conjecture remains elu- 
sive, all existing evidence points at some version of it be- 
ing true in generic 4-dimensional, asymptotically flat space- 
times 101. This is reassuiing because it means that a quan- 
tum theory of gravity is not needed to understand the astro- 
physics of gravitationally-collapsed objects^, since all curva- 
ture singulaiities would be cloaked behind an event horizon 
and would therefore be inaccessible to external observers. 

Because, according to the Cosmic Censorship Conjecture, 
curvature singularities are genetically hidden behind black- 
hole event horizons, the most natural way to search for naked 
singularities is to look for instabilities of the Kerr geometry, 
which is the most general stationary vacuum black-hole solu- 
tion of Einstein's equations in a four-dimensional, asymptot- 
ically flat spacetime |9(1. At the linearized level, all existing 
evidence points to the stability of the Kerr event horizon'': 

(i) the Kerr metric outside the event horizon is perturba- 
tively stable against exponentially growing modes of massless 
fields iflTll . The Kerr exterior geometry is unstable against 
massive scalar fields lfl2ljl7ll . but the instability is thought to 
extract rotational energy from the geometry, and does not in 
principle destroy the horizon. 

(ii) Because naked singularities appear in the Kerr geometry 
when the angular momentum of the spacetime is above the 
Kerr bound, i.e. for cJ /GM^ > 1, a simple gedanken ex- 
periment was proposed by Wald ifisll to try to create naked 
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singularities. This consisted of throwing point particles hav- 
ing large angular momentum into an extremal Kerr black hole. 
In the test-particle approximation {i.e. at the lowest order in 
the mass of the particle), Wald finds that particles that would 
over-spin the geometry past the Kerr bound are not captured, 
but are scattered. This conclusion generalizes to a wide vari- 
ety of backgrounds lfl9ll . 

On the other hand, Kerr spacetimes that present 
cJ/GAP > 1, and which therefore contain naked singular- 
ities, are linearly unstable ||20|| . Thus, a consistent picture 
emerges, suggesting that black holes are stable and do not give 
rise to naked singularities. 

At the full nonlinear level, all numerical-relativity evolu- 
tions of Einstein equations that asymptote to the Kerr geom- 
etry at late times are consistent with this scenario, suggesting 
stability of the Kerr spacetime ||2T1 - I28ll . In particular, a pos- 
sible generalization of Wald's process to comparable mass ra- 
tios has been studied in Refs. ll2lll22ll . which considered high- 
energy collisions between two comparable-mass non-spinning 
black holes. Refs. ll2lll22ll found no evidence of formation of 
naked singularities, because either the black holes are simply 
scattered to infinity, or the full nonlinear equations make the 
system radiate enough angular momentum to form a single 
black hole. (More specifically, for finely tuned impact param- 
eters, the black holes spend a finite amount of time in zoom- 
whirl orbits and radiate the excess angular momentum, even- 
tually being "allowed" to merge.) 

Similarly, Abrahams et al. ll29tl studied the collapse of ax- 
isymmetric tori made of collisionless matter, and found that 
tori with cJ / GAP < 1 ("sub-Kerr" configurations) collapse 
to Kerr black holes, while "supra-Kerr" tori with cJ /GAP > 
1 collapse to new equilibrium configurations. The collapse 
of rotating stars has also been investigated by a number of 
authors |[30l - [35ll . These studies show that the collapse of sub- 
Kerr stellar models produces Kerr black holes, while the col- 
lapse of supra-Kerr models does not give rise to naked sin- 
gularities. For instance, Ref. jssil shows that the collapse of 
a supra-Kerr differentially rotating polytropic star produces a 
torus, which then fragments in clumps that merge again, form- 
ing a bar and eventually a stable axisymmetric configuration. 

Recently however, Jacobson and Sotiriou (JS) |l36|| (build- 
ing on Refs. IstII ) have re-considered Wald's process, but us- 
ing almost extremal black holes, rather than exactly extremal 
ones. Surprisingly, JS showed that this change is enough to 
allow test-particles with dangerously large angular momen- 
tum to be captured, over-spinning the black hole and creating 
naked singularities. However, the particles considered by JS 
need to have a finite mass-energy to overspin the black hole, 
and in this sense the situation that they consider is interme- 
diate between Wald's original construction for test particles, 
and the full nonlinear analysis of Refs. ll2lll22ll for compara- 
ble masses. Therefore, as acknowledged by JS themselves, a 
test-particle analysis such as that of Wald's gedanken exper- 
iment is likely not adequate for this scenario, because it ne- 
glects the conservative and dissipative self-force, which may 
be important Issll . 

Here we expand on our previous Letter js^, and show 
that the dissipative self-force (equivalent to radiation reaction. 



i.e. the energy and angular momentum losses through gravita- 
tional waves) can prevent the formation of naked singularities 
only for some of JS's orbits. However, we will show that for 
all these orbits the conservative self-force is comparable to the 
terms giving rise to naked singularities, and should therefore 
be taken into account. 

In particular, in Sec. we review the orbits that JS iden- 
tified as giving rise to naked singularities. In Sec. |III] we 
show how to calculate the gravitational-wave fluxes produced 
by these orbits, both analytically (Sec. lIII AT ) and numerically 
(Sees. HITB] and HITCl ). A comparison between the analytical 
and numerical fluxes is presented in Sec. |IV] while in Sec. |V] 
we analyze the effect of the conservative self-force on the JS 
process. In Sec. |VT] we present our conclusions. Hereafter we 
set G = c = i\/ = 1. 



II. THE OVERSPINNING ORBITS OF JACOBSON & 
SOTIRIOU 

Let us consider a Kerr black hole with spin parameter 
a = J /AP = 1 — 2e^, with e ^ 1, and a non-spinning 
test-particle with energy E, angular momentum L and mass 
/i. By adopting suitable length units, we can set the mass AI 
of the black hole to 1 without loss of generality, and it there- 
fore follows that in order for the test-particle approximation 
to be valid it must be /^j <C 1 as well as <C 1 and L ^ 1. 

Neglecting the loss of energy and angular momentum 
through gravitational waves (which is equivalent to the so- 
called dissipative self-force) and the conservative self-force 
{i.e., the modifications of the effective potential for the parti- 
cle's motion due to the small but finite mass ratio of the sys- 
tem), the particle moves on a geodesic of the background Kerr 
spacetime. JS then impose f /) that this geodesic orbit fall into 
the black hole, and f /;) that the Kerr black hole be spun up past 
the extremal limit when the particle is captured, i.e. that the 
final spin aj? ^ — {a + L) / {1 + E)"^ be larger than 1. Condition 
( i) implies an upper limit L < Lmax on the angular momen- 
tum, because particles with large angular momenta are simply 
scattered, while condition ( ii) implies a lower limit L > Lmin 
on the angular momentum, because the particle needs to trans- 
fer a sufficient amount of angular momentum to the black hole 
in order to spin it up past the extremal limit. In particular, one 
finds 

L„,in = 2e^ + 2E + E^ <L< L,nax = (2 + 4e)E . (1) 

In order for orbits satisfying both ( i) and ( ii) to exist, one has 
then to impose Lmax > ^min, which yields 

Enun - (2 - V2)e <E< E„,,^ = (2 + V2)e . (2) 

This interval is indeed finite when e > 0, while it shrinks to 
zero when e = {i.e., when a = 1). This confirms Wald's 
results that an extremal black hole cannot be spun up past the 
extremal limit by the capture of test particles. Lastly, JS no- 
ticed that intervals ^ and (|2]i contain both bound orbits {i.e. 
orbits that start with zero radial velocity at finite radius) and 
unbound oihits {i.e. orbits that start from infinity). 
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Parameterizing the above intervals as 

E = E'min + ^(E'niax " Ea^in) = £^min + 2x^26 (3) 
L = Lmin + y{Lmax ^ imin) = ^min + 8ye^(l - x)x (4) 



with 0<a:;<l,0<y<l, the final spin is 

a + L 



Uf — 



{1 + Ey 



= l+8e^(l-a;)xy + 0(e-^)>l, (5) 



which is indeed larger than 1 by terms that are quadratic in 
e. It is clear, however, that the inclusion of the gravitational- 
wave losses of energy and angular momentum, Siad and irad. 
can in principle affect JS's analysis by changing the prediction 
(|5]l for the final spin to 



a/ = 1 + 86^(1 - x)xy + 2E, 



ad 



-'rad 



+ O(e') . (6) 



In order to analyze the magnitudes of Siad and iiad relative 
to the spin-up term in (|6]l, we restrict our attention on unbound 
orbits,'* and following JS we assume E/i^i ^ 1 and ^ 1 
(null orbits) to further simplify our analysis. These relativistic 
orbits are characterized by the impact parameter b = L/E 
alone. From Eqs. ([TJ and it follows that JS's orbits have 
L ^ bE, withfe = 2 + 4e[l-2x{x-l){y-l)/{2 + V2{2x- 
1 ) ) ] . Varying x and y between and 1 , one obtains b = 2 + Se 
with 2\/2 < (5 < 4. However, because 6ph = 2 + 2\/3e + 
O (e^ ) is the impact parameter of the circular photon orbit (or 
"light-ring"), only orbits with 2%/2 are unbound. 

When S w 2-^3, the particle has an impact parameter very 
close to that of the light-ring, and therefore orbits around the 
light-ring many times before plunging into the black hole. The 
emission of gravitational waves should be important for these 
orbits (because the gravitational-wave fluxes will be propor- 
tional to the number of cycles at the light-ring), and could 
prevent the formation of naked singularities or at least invali- 
date JS's analysis. In fact, for S arbitrarily close to 2\/3, the 
particles would orbit around the light-ring an arbitrarily large 
number of times, and gravitational- wave emission must be im- 
portant |40l]. We will show, however, that this is not true for 
all of JS's orbits but only for a subset of them. Nevertheless, 
we will also show that for all of JS's orbits, the conservative 
self-force is always important and seems to have the right sign 
to prevent the formation of naked singularities. 

To estimate the magnitudes of E'rad and iiad, we first com- 
pute the number of cycles A'cycios described by the relativistic 
particle at the light-ring using the geodesies equations. As 
we have mentioned, E'lad and Liad are expected to be propor- 
tional to iVcycics, because the plunge from infinity to the light 
ring and that from the light-ring to the horizon are not ex- 
pected to produce significant amounts of gravitational waves, 
since they happen on a dynamical timescale. 



As we mentioned, JS also considered bound orbits, falling into the Kerr 
black hole from a Boyer-Lindquist radius r = rj^o^ + 0{e) (rjjo,. being the 
horizon's radius). However, these orbits pose a problem, as we will show 
later, because the distance to the horizon is comparable to the particle's 
minimum attainable size max(i?, /^) > e, so finite-size effects should be 
taken into account. 



The geodesic equations for null equatorial orbits read 

{dr/d\f ^Vr{r), d(j)/d\ = V^{r), (7) 
Vr{r) =r[r^ - Aah - b^{r - 2) + a^{r + 2)] , (8) 
[2a + b{r^2)]r 



a2 + (r-2)r ' 



where A is a non-affine parameter (dX = dXamncE/r'^). For 
b = &ph(l — k), with fc ^ e <C 1, the radial potential presents 
a minimum at r = r„dn: 

T/'"'"«8fce/V3 + 0(fce2) + 0(fc2), (10) 
r-min ~ '^ph - 4fc(l + V3e)/3 + 0{ke^) + 0{k^) . (11) 

Near this minimum one has 

Vr{r) = V;,"""+ iK"(^mi„)(r - r„dn)'+ 0{r - r,„i„)3 (12) 
V^r) ^^ + ^+0{e)+0{k)+0{r- r,„i„) , (13) 
where V^'{rnun) = 6 + O (e) + O (fc). Therefore, we get 

-1/2 




—^ke + 3(r - n-ninf 



(14) 



Integrating from rmin — Ari to rmi,i + Ar2, with Ari 2 S> ke, 
the number of cycles near the minimum is 



N, 



cycles 



i+Ar2 



,-Ari 



d4> dr 
dr 27r 



[A+B log (ke)] 



(15) 

A and B being constants depending on the integration inter- 
val. Fixing e, and thus the black-hole spin, we can see that 
^cycles depends on log k, and diverges when k ^ ll40ll . 

In Sec. HV] we will study gravitational-wave emission nu- 
merically for geodesies having E = (i?max + £'min)/2 = 2e, 
L = bp\^E{l — k) with k = ji = 10^^, in Kerr spacetimes 
with a = 0.99, 0.994, 0.998 and 0.9998. In particular, we in- 
tegrate these geodesies numerically (using a 4'^ order Runge- 
Kutta integrator with Richardson extrapolation). The Boyer- 
Lindquist radius a function of coordinate time is shown in 
Fig- [11 As can be seen, after an initial plunge, the orbits show 
a plateau (corresponding to the particle orbiting at the light- 
ring) followed by the final plunge towards the horizon. The 
time spent at the light-ring increases with a {cf. Eq. (fTStl. 
and the horizon radius is approached in an infinite coordinate 
time, in agreement with the analytical behavior pointed out 
by Ref. ll41fl . In order to estimate A^cycios, we count the num- 
ber of cycles of these numerically-integrated geodesies from 
r = 1.05rph to r = (1.9rph + 0.1rhor)/2 (rph being the 
light-ring radius). These values for A^cycies can be obtained 
from Eq. O with A « 0.620, B « 0.01515, to within 2.5% 
(- 0.17 cycles). 



4 




a=0.99 
a=0.994 
a=0.998 
a=0.9998 



200 400 600 800 1000 1200 1400 
t 

FIG. 1 : The radial evolution with coordinate time, for geodesies hav- 
ing E = (S,„ax + Emin)/2 = 2e and L = 6ph-B(l - k) with 
fc = = 10"^, in Kerr spacetimes with a = 0.99, 0.994, 0.998 and 
0.9998. In particular, r and rn are the Boyer-Lindquist coordinates 
of the particle and of the horizon. 



III. THE GRAVITATIONAL-WAVE FLUXES 
A. Analytics 

The radiation emitted by point particles on circular 
geodesies in the Kerr geometry was studied semi-analytically 
by Chrzanowski and by Chrzanowski & Misner, among oth- 
ers ll42l - l45ll . The original motivation for these investigations 
was the possible existence of synchrotron radiation in the 
Kerr geometry, a possibility that was ruled out by these stud- 
ies. Performing a WKB analysis of the Teukolsky equation, 
Chrzanowski ll42ll and Chrzanowski & Misner llisll concluded 
that the energy flux by a particle on the innermost stable cir- 
cular orbit (ISCO) is suppressed as the geometry approaches 
the extremal Kerr geometry. In particular, they found that the 
flux and therefore the total energy radiated by a particle with 
mass p, scales as 



{r - rH)^-l^ 



(16) 



where the orbital radius r approaches the horizon radius tr 
in the extremal limit (because risco = fH for a = 1, cf. 
Ref. S). 

This suppression of the energy flux in the extremal Kerr 
limit is common to other massless fields, like scalar or elec- 
tromagnetic fields 1I42I - I44I1 . A possible interpretation for this 
behavior was put forward by Breuer ||44| . who argues that in 
the a 1 limit, a particle close to the ISCO approaches a 
principal null direction of the Kerr geometry, for which tidal 
and radiation effects are strongly suppressed lliTllisl . 

Alternatively, Eq. ( fT6] l can be derived from dimensional ar- 
guments, which suggest that dE/dr ^ p^ [r being the parti- 
cle's proper time; see also Eq. dTTTil. with the factor {r — rn) 
coming from the gravitational redshift at the ISCO. In fact, 
for almost extremal black holes the ISCO is very close to 
the horizon, and for near-horizon orbits one has dt/dr ~ 



TH/{f — "Th) 114 ill . This second derivation suggests that the 
scaling ( fT6] l should be valid also for the fluxes falling into 
the black-hole horizon, and not only for the fluxes at infinity 
studied by Refs. iH-Q. In fact, Kesden (H] has recently 
revisited the problem of studying the gravitational emission 
by a particle on the ISCO of a Kerr spacetime. He used the 
frequency-domain Teukolsky code GREMLIN ll49ll . which 
can calculate both the fluxes at infinity and those falling into 
the horizon, and confirmed that the scaling (fTSI l holds not only 
for the outgoing fluxes but also for the total ones. 

Because the ISCO approaches the light ring in an almost 
extremal Kerr spacetime (the ISCO, light-ring and horizon ac- 
tually coincide in the extremal case ll46ll ). it should be possible 
to use the Eq. ( fT6b to infer the scaling of the gravitational radi- 
ation emitted by a photon at the light-ring. The same scaling 
will then hold also for ultrarelativistic particles (i.e. ones with 
E/ p ^ 1) that orbit near the light-ring many times, such as 
those in which we are interested in this paper (cf. Sec |II}. 
This is because any null geodesic (in a generic geometry) can 
be obtained as the limit of a series of timelike geodesies for 
which the particle's 4-momentum is kept finite and the parti- 
cle's rest-mass goes to zero. [More intuitively, this amounts 
to saying that in General Relativity the (positive-energy) or- 
bits of a neutrino, which has a small but non-zero mass, are 
very similar to those of a photon.] 

Therefore, in order to take the null-geodesic limit of 
Eq. ( fTSI ), we first have to re-express it in terms of the energy 
E of the (massive) particle at the ISCO. Because for such a 
particle E ^ p,we can recast Eq. ( fT6] l in the form 



Eir^ir- rH)E^ 



(17) 



Because when a ^ 1 the ISCO radius approaches the light- 
ring, this scaling must be also valid for a particle on an un- 
stable circular orbit (such orbits exist between the ISCO and 
the light-ring ifsoll ). The energy of the unstable circular orbits 
becomes positive near the light-ring (more specifically, it is 
positive for orbits with radius Tph < r < r^^, where is 
the radius of the marginally bound circular orbit (Bot ). We can 
then take a series of unstable circular orbits with constant pos- 
itive energy E and with radius approaching the light-ring. The 
rest-mass of these orbits will clearly go to zero as r — > rph, 
because the specific energy E/ jjL must diverge at the light- 
ring 1I50I1 . Therefore, this series of timelike geodesies satisfies 
the conditions, outlined above, necessary to approach a null 
geodesies. This null geodesies will of course be a circular or- 
bit at the light-ring with the same energy E as the sequence 
of unstable circular orbits that we have just mentioned, and its 
flux must too be desciibed by the scaling ( fTTj i. 

We stress that this result makes sense physically because 
in General Relativity it is the energy and not the rest-mass 
that gravitates. In fact, a more direct deiivation of Eq. dTT] ) 
for a photon at the light-ring should be possible by solving 
the Teukolsky equations using the stress-energy tensor for a 
photon, which is 



P^(AK(A) 



(5W(x-z(A)) 



(18) 
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where z(A) is the wordUne and p'^ = dz^/dX is the 4- 
momentum. (See for instance Ref. llsill for a derivation of 
the stress-energy tensor of the electromagnetic field in the 
geometric-optics limit.) Clearly, T^^ is non-zero despite the 
mass of the photon being zero, because the photon does have 
an energy, which curves the geometry. Of course, T^^ agrees 
with the stress-energy tensor of massive particle [cf. Eq. 
(19.3) of Ref. [52]], when the mass is sent to zero keeping 
the 4-momentum finite. 

Assuming then the validity of Eq. ( fTTj i for ultrarelativistic 
orbits and for photons as a working hypothesis, to be con- 
firmed numerically later in this paper, we can explore its con- 
sequences for the JS orbits considered in the previous section. 
For those orbits, E ^ r — rn ^ £ and we therefore expect 



(19) 



Also, we can now derive the behavior of the total radiated en- 
ergy for high-energy plunges with near-critical impact param- 
eter Because for these orbits the particle spends most of its 
time on almost circular orbits at the light ring, one can write 

Erad = AE{e) X A'^cyclos , irad = AL(e) X iVcyclos , (20) 

where AE and AL are the fluxes in a single orbit. From 
a frequency-domain analysis, AE/AL must equal the hght- 
ring frequency, fiph w 1/2 — (\/3/2)e, hence 

AE{e) = £;i(e)(l + e2e), (21) 
AL(e) = 2£;i(e)[l + (x/3 + e2)e]. (22) 

Here 62 is an undetermined coefficient. 

B. Numerics 

The Teukolsky master equation describes scalar, vector and 
tensor field perturbations in the space-time of Kerr black 
holes ifssll . In Boyer-Lindquist coordinates, this equation 
takes the form 



Jl\2 



c? sin^ I 



-2s 

-A-' 
1 



A 

dr (A'^+^a.^-; 
2- 



ia cos 9 
1 



sm C/ L 

- (s^ cot^ e 



2s 



-deiainedg^) 
u 

a{r — 1) ico 



-47r(r^ 



cos 



2q\ 



(23) 



where A ~ — 2r + a^ and s is the "spin weight" of the field. 
The s = — 2 version of this equation describes the evolution 
of the Weyl scalar -04 = ^'/(r — ia cos OY that characterizes 
the outgoing gravitational radiation. This is sufficient for the 
present work, because in order to determine the final spin a j 
one only needs to account for the gravitational-wave radiation 
that leaves the binary system, while the fluxes emitted by the 



particle into the black-hole horizon ("ingoing fluxes") cancel 
out when computing the total angular momentum and energy 
and therefore do not affect the final spin parameter^ 

Computing the radiative energy and angular momentum 
loss is fairly straightforward from 'i/. More specifically, we 
use the following expressions for the radiated energy and an- 
gular momentum fluxes 



dE 
'dt 

dt 



= lim 



1 



dn 



= — lim 



Re 



t 2> 

dt>(i>,^,(^) 



dn a 



r— >oo 

XI/ dt' I di'^{i,r,e,ip) 




The variable T on the right hand side of Eq. ( |23] | is a point- 
particle source-term constructed from the energy-momentum 
tensor describing a particle moving in a Kerr spacetime. The 
particle energy-momentum tensor has the form 



Tap 



Etsine 



5[r^r{t)]5[e-e{t)]6[^-cj>{tm^) 



where is the 4-velocity of the particle. It is noteworthy 
that t = dt/dr appears in the denominator of this expres- 
sion. As the particle approaches the horizon, i — > 00, which 
is just the well-known "infinite redshift" effect at the horizon 
of a black hole. Thus the particle source-term smoothly "red- 
shifts away" as the particle approaches the horizon, therefore 
allowing the Teukolsky equation to gradually transition into 
its homogeneous form. This smoothly and naturally connects 
the gravitational radiation from the last few orbital cycles to 
the black hole's quasi-normal modes. 

To solve Eq. ( |23] ) numerically in time-domain we take the 
approach first introduced by Krivan et al. in Ref. t54j] . Our 
code that solves the Teukolsky equation uses the same ap- 
proach, therefore the contents of this section are largely a re- 
view of the work presented in the relevant literature. 

Our time-domain code uses the tortoise coordinate r* in 
the radial direction and azimuthal coordinate 0. These coor- 
dinates are much better suited for performing numerical evo- 
lutions, as detailed in Ref. [54]. They are related to the usual 
Boyer-Lindquist coordinates through the equations. 



dr* 



A 



-dr 



and 



-dr 



(26) 



(27) 



Following Ref. 115411 . we factor out the azimuthal dependence 
by using the m-mode decomposition. 



e™'^r3$(t,r*,6'). 



(28) 



^ It is conceivably possible, however, that the ingoing fluxes might overspin 
the black hole and possibly create a naked singularity before the particle is 
captured. We will comment on this effect in Sec. lIVI 
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Defining 

n EE dt<i> + bdr*^ , (29) 
6^^, (30) 

wliere 

= {r^ + - A sin^ e (31) 

allows the Teukolsky equation to be rewritten in first-order 
form 

dtu + Mdr*u + Lu + Au = T, (32) 

where 

M = {$fi,*/,nfi,n/} (33) 

is the solution vector The subscripts R and / refer to the 
real and imaginary parts respectively (recall that the Teukol- 
sky function ^l* is a complex valued quantity). Explicit forms 
for the matrices M, A and L can be easily found in the rele- 
vant literature ll54ll . 

Lastly, an explicit time-evolution numerical scheme is de- 
veloped for this first-order, linear PDE system using the two- 
step, 2nd-order Lax-Wendroff finite-difference method. Ex- 
plicit details on this approach can be found in Ref. llsill . Sym- 
metries of the spheroidal harmonics are used to determine the 
angular boundary conditions: For even |m| modes, we have 
de^ = at 6* = 0, TT while $ = at = 0, tt for modes of 
odd |m| . We set $ and 11 to zero on the inner and outer radial 
boundaries. 

One major challenge in numerically solving Eq. (|23] | is 
representing a point-like particle source on a numerical grid. 
There are multiple approaches towards tackling this problem, 
such as representing the particle as a narrow Gaussian distri- 
bution IdSII or taking a more efficient "discrete Dirac-delta" 
approach as presented in Ref. [:56!|. One observation that we 
make specific to this work is that for near-extremal black holes 
in Boyer-Lindquist coordinates, the particle's orbit, the light- 
ring and the horizon are extremely close, therefore we mod- 
eled the point-particle to have a fixed width in the tortoise 
coordinate r* as opposed to r. This allows our code to resolve 
these distinct regions, because they are relatively widely sep- 
arated in r*. To test that our entire numerical scheme is work- 
ing properly we performed extensive convergence tests on the 
data generated by our code, especially for the near-extremal 
Kerr hole cases. For these convergence tests, we kept the ratio 
of the particle's width to grid spacing constant i.e. upon dou- 
bling grid-density we consistently halved the particle width. 
Sample convergence results for the a — 0.9998 and mode 
771 — 2 case are presented in Fig. |2] We observe a conver- 
gence rate that is extremely close to the expected 2nd-order 
convergence. 

One possible complication in our analysis comes from the 
fact that the convergence rate is actually slightly faster than 
the theoretical 2nd-order convergence, as can be seen from 
Fig. [21 This makes a Richardson extrapolation to further re- 
duce discretization errors difficult to implement (because one 
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FIG. 2: Demonstration of our time-domain Teukolsky code's conver- 
gence rate, for a = 0.9998 and tu — 2. 



should account for the changes of the convergence rate along 
the trajectory in order to obtain an accurate result). There- 
fore, we simply decide to decrease the grid spacing and par- 
ticle's width until our numerical results change by no more 
than 0.5%. The final grid spacing that we use in this work 
is 1/100 in the r*-coordinate and 0.02 radians in the ^-grid, 
and we use 1 /500 as the time-step. We also verify a posteriori 
that our numerical errors are less than 0.5% by checking that 
the ratio AE/AL coincides with the light-ring frequency to 
within 0.35%. 

Our code that implements this numerical scheme is a 
Fortran-code that is fully parallelized to execute efficiently 
on a computer cluster The parallelization approach taken is 
the standard domain-decomposition approach (on the radial- 
coordinate numerical grid) with message-passing enabled us- 
ing OpenMPI. Good scaling has been observed for several 
hundred processor-cores. In this current work, we made use 
of 200 processor-cores for computing each m-mode for every 
case that we studied. 



C. Summing over m-modes 

As noted above, our evolution code evolves each m-mode 
separately. Therefore, we obtain the fluxes and radiated en- 
ergy and angular momentum for each 77i-mode through dis- 
tinct and separate numerical simulations. In order to com- 
pute the total radiated quantities, we would need to perform a 
sum over the results obtained from these different simulations. 
However, for higher m-modes it becomes increasingly diffi- 
cult to perform accurate numerical computations mainly due 
to the requirement of significantly higher resolution (higher 
?T!-modes involve higher values of i that require higher angu- 
lar grid resolution for accurate representation; there is also a 
Courant condition that requires us to reduce the time-step with 
higher angular resolution, thus making the computation even 
more demanding). Therefore, we use an alternative approach 
to estimate the radiated quantities from m-modes higher than 
m — 10. For the case of circular and equatorial orbits, Finn & 
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Thome ll57ll show that 



Em — 



2(m + l)(m + 2)(2m + 1)! m' 
(m- l)[2"m! (2to+ 1)!!]2 



2m+l 



(34) 

where i7 in Eq. ( [34l i is the Keplerian angular frequency of the 
orbit i.e. fi^/^ = 1/r. 

Because the gravitational-wave fluxes for JS orbits are pro- 
duced, for the most part, during the quasi-circular equato- 
rial cycles at the light-ring, our numerically calculated fluxes 
should satisfy this scaling, and in particular they should be in 
a geometric progression for large values of m, because Eq.[34l 
implies Em+i/Em — > constant for m ^ 1. We have in- 
deed verified that this is the case, i.e. we have checked that 
Em.+i/Em and L„i+i/i,„ are constant to within 0.5% al- 
ready for m > 9. 

This behavior allows us to reconstruct the total fluxes in the 
following way: (1) We use our Teukolsky evolution code de- 
scribed above to compute the fluxes and the radiated quantities 
for TO-modes up to m = 10; (2) We use the fluxes calculated 
for m = 9 and rn = 10 to estimate the ratio entering the 
geometric progression described above and thus obtain an es- 
timate of the contribution from the higher m-modes; (3) We 
finally add in this estimate to the explicitly computed values 
in step (1) and thus obtain the total radiated amount for both 
the energy and the angular momentum. 

We stress that while the fluxes summed up to m = 10 are 
accurate at least to within 0.5% (cf. discussion in the previ- 
ous section), the procedure that we have just described intro- 
duces larger errors into the total fluxes. For the orbits that 
we consider in this paper, the asymptotic ratios Em.+i/Em 
and Lm+i/Lm grow from 0.928 for a = 0.99 to 0.94 for 
a = 0.9998. Since the sum of a geometric series is propor- 
tional to 1/(1 — r), r being the asymptotic ratio, a 0.5% er- 
ror in estimating r would lead to a ~8% error in the sum of 
the fluxes with m > 10. This propagates into a significant 
error in the total fluxes, because for relativistic plunging or- 
bits around almost-extremal Kerr black holes such as those 
that we are considering in this paper, the contribution of the 
large multipole moments to the total fluxes decays slowly with 
TO. In particular, because the contribution of the fluxes with 
TO > 10 grows from 1.25 times the contribution of the fluxes 
with m = - 10 for a = 0.99, up to 1.7 times for a = 0.9998, 
we estimate the error of our total fluxes to be < 5%. 



COMPARING THE ANALYTICAL AND NUMERICAL 
GRAVITATIONAL-WAVE FLUXES 



IV. 



In our previous Letter 13911 . we considered black holes with 
a = 0.99, 0.992, 0.994, 0.996, 0.998, 0.999 and 0.9998, 
and geodesies having E = (i?max + £'niin)/2 = 2e, L = 
6ph£'(l - k) with k = 10"^ and ^ = 0.001. Using the 
time-domain Teukolsky code that we described in Sec. IIIIBI 
we verified that the theoretical scaling given by Eqs. (fT9t-(l22]i 
works well for a < 0.999, but predicts fluxes that are too 
large for a = 0.999 and a = 0.9998. More specifically, for 
a < 0.999 the deviations of the numerical results from the 



theoretical scaling are about 1 — 3% and therefore comparable 
to the errors discussed in Sec. IIIIBl However, for a = 0.999 
and a = 0.9998 the fluxes predicted by Eqs. ([T9]l-(|22ll are 
respectively 12% and 84% larger than the numerical ones. 

Because our Letter li3^ was only concerned about whether 
the gravitational-wave fluxes could prevent naked singulari- 
ties from forming, this discrepancy did not affect our conclu- 
sions. In particular, it reinforced our finding that there are 
orbits giving rise to naked singularities even when radiation 
reaction is taken into account. In this section, however, we 
will investigate the origin of this discrepancy between our the- 
oretical scaling and our numerical results. 

First, let us note that a crucial assumption in the deriva- 
tion of the scaling (fT9]l-(l22]i is that the orbits be relativistic. 
In fact, only under that assumption we can expect the fluxes 
to depend only on the energy E (and therefore on e) and not 
on the rest mass fi. However, because the orbits considered 
in Ref. 13§] have energy E = (-Emax + -Ei„in)/2 = 2e and 
mass fi = 10~3, the ratio E/^i is - 22 for a = 0.999 and 
10 for a = 0.9998. Therefore, a possible explanation for the 
disagreement between our analytical and numerical results is 
that the orbits that we considered are simply not relativistic 
enough. In order to test this hypothesis, we have tried produc- 
ing orbits that have exactly the same energy and angular mo- 
mentum as those considered in Ref. [39], but which have dif- 
ferent rest mass. More specifically, we have considered black 
holes with a = 0.99, 0.992, 0.994, 0.996, 0.998, 0.999 and 
0.9998, and geodesies having E = (E'max + ^^min) /2 = 2e, 
L = feph£'(l - k) with k = fj. = 10"^. If the orbits with 
ji = 10^""* were already sufficiently relativistic to satisfy the 
scaling (fT9]l-(l22li. then their fluxes should coincide, to within 
the numerical errors, with those of the orbits with /^t = 10^^, 
since the two set of orbits have exactly the same energies and 
angular momenta. 

In Fig. [3j we therefore show the fractional difference be- 
tween the energy fluxes (in a single orbit at the light ring) for 
the two sets of orbits. In particular, we have calculated this 
difference both for the total fluxes (obtained by summing over 
all multipole moments to as described in Sec. IlllCt and for 
the "partial" fluxes obtained by summing up to m = 10. As 
can be seen the difference is less than 0.5% for a < 0.999, is 
about 1% for a = 0.999 and grows to -10% for a = 0.9998, 
thus suggesting that at least in this last case the orbit with 
/i = 10^"^ is not sufficiently relativistic. As a further con- 
firmation, we have produced numerical fluxes for a particle 
moving in a Kerr spacetime with a = 0.9998, and having the 
same energy and angular momentum as the previous orbits but 
rest mass fi = 10^^. As can be seen from Fig. |3] this orbit 
gives fluxes that agree with those produced for the /i = 10~^ 
orbit to within less than 0.1%. 

This test confirms that orbits with /i = 10~^ are sufficiently 
relativistic for the scaling (fT9t-(l22t to hold. However, al- 
though using these orbits improves the agreement with the 
expected scaling, the deviations remain as large as 5% for 
a = 0.999 and 70% for a = 0.9998. In fact, a careful analysis 
of the convergence of all our numerical fluxes with respect to 
the grid spacing and particle's width revealed that the modes 
with TO > 8 did not completely converge in the case of the 
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orbits around black holes with a = 0.999 and a = 0.9998 
considered in our previous work ll39l . This convergence error 
was then amplified by the procedure that we use to reconstruct 
the contribution of the large-m modes. (As we discussed in 
Sec. unci this procedure is very sensitive to the numerical 
results for m = 9 and m = 10). 

The resolution that we used was instead sufficient in the 
case of the other multipole moments for a = 0.999 and 
a ~ 0.9998, and for all the multipole moments at lower spins. 
Indeed, it is not surprising that this problem affected only the 
high-m modes for spins very close to 1. When a ^ I, the 
frequency of the light-ring becomes very close to that of the 
horizon, and so does the radius (in Boyer-Lindquist coordi- 
nates). It is therefore necessary, as already mentioned in Sec. 
IIIIBI to use a very high resolution in the near-horizon region 
to calculate the fluxes accurately for m ^ 1 (because large 
to's correspond to small lengthscales). 

In order to amend this problem, we focused on just four or- 
bits - namely orbits around black holes with a = 0.99, 0.994, 
0.998 and 0.9998, and having E = (£;,„ax + £^min)/2 = 2e, 
L = bpiiE{l — k) with k ~ — 10~^ - but calculated 
the fluxes with very high grid resolution and small particle's 
width. More specifically, as explained in Sec. IIIIBI we es- 
timated the error connected to the finite grid resolution and 
particle's width to be less than 0.5%, and the error due to the 
reconstruction of the large-m modes to be less than 5%. 

As discussed in Sec. [Ill the number of cycles described by 
these orbits at the light ring is reproduced by Eq. (ITSt with 
A « 0.620 and B w 0.01515 to within 2.5%. Assuming 
then El = eie", we fit the numerical energy and angular- 
momentum fluxes with Eqs. (|20]|-(|22]|. obtaining n w 2.95. 
Because this value is very close to the theoretical value n 3 
{cf. Eq. (fT9l)1. we then assume n ~ 'i and fit the data with 
only two free parameters, ei and 62, obtaining ei = 233.72 
and 62 = —5.83.^ With these values, Eqs. (fT9]l-(l22]i repro- 
duce the numerical data with residuals < 4%. As discussed 
in Sees, nil Bj and HITCl these residuals are comparable to the 
errors affecting our total fluxes, thus confirming our expected 
scaling. 

Utilizing Eq. it is now straightforward to determine 
the consequences of this scaling for the final spin a/. From 
Eqs. (I20] i-(l22li. we obtain 

iiad - 2£'i.ad = 2\/3e£;iA'^cyclcs , (35) 

which implies in particular Lj-ad — 2iJi.ad > 0, i.e. the 
fluxes tend to decrease the final spin. Moreover, using expres- 
sion (flSl l for A'cycios ™d expression ( fT9] l for Ei, we obtain 



* This value of ei, which represents the overall normalization of the fluxes 
emitted in a single orbit, is significantly different from that reported in 
Ref. (3^. This is due to the different ways of calculating the number of 
orbits jVpycicB at the light ring that we used in this paper and in Ref. I39I1 . 
In Ref. 13911 . we defined Afcydcs the number of cycles between r = 
l.OSrpij and r = {r^i^ + rijoj)/2, while here we decided to bracket the 
light ring more closely and considered the number of cycles between r = 
l.OSrpij and r = (1.9rpij + 0.1rijoj)/2. These different definitions are 
also the reason why the values of A and B that we report in this paper are 
significantly different from those of Ref. js^. 



|Ae(h=io-=)/Ae(h=io-^)-i| (total) — 




0.0001 I ' ' ' ' 1 
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a 

FIG. 3: Fractional difference between the energy fluxes emitted in a 
single orbit at the light ring for fj, = 10""^ and fi — W^^ (filled sym- 
bols), and for fi = 10^^ and fi — W~"' (empty symbols). The fluxes 
are obtained by summing over all multipole moments ("total" fluxes, 
represented by squares) using the prescription outlined in Sec. lIIICl 
or by summing only up to m = 10 (fluxes "up to m = 10", repre- 
sented by circles). 



a 


0.99 


0.994 


0.998 


0.9998 


aj 


0.886 


0.963 


0.997 


1.00004 




1.0043 


1.0026 


1.0009 


1.00009 



TABLE I: Initial and final black-hole spin after absorbing a relativis- 
tic particle with energy E — \/2{l — a) and angular momentum 
L = feph-E(l — 10~^), neglecting conservative self-force effects, but 
not radiation reaction. We also show the final spin without radiation 
reaction (aj^) predicted by JS. 



that Lrad ~ 2i?i.ad ^ e'' log(fce). Comparing this scaling with 
the "overspinning" term 8e^(l — x)xy in Eq. (|6]l, it is clear 
that if k is kept constant and e is sent to {i.e., if one con- 
siders initial spins larger than some critical value Ocrit), the 
effect of the fluxes eventually becomes subdominant relative 
to the "overspinning" term. In other words, the fluxes tend to 
decrease the final spin a/, but their effect is not sufficient to 
prevent the formation of naked singularities if the initial spin 
is sufficiently close to 1 . Our numerical results confirm this 
picture, as can be seen from Table |I] where we show the ini- 
tial spin a, the final spin aj^ as computed by JS (i.e., without 
including the effect of the fluxes), and the final spin af, tak- 
ing into account radiation reaction. As can be seen, for the 
orbits that we consider aj > 1 already for a = 0.9998, and 
therefore the critical spin mentioned above is Ocrit ~ 0.9998 
(corresponding to e w 0.01 |[39ll ). 

Finally, let us further comment on the effect of the ingoing 
fluxes. As already mentioned in Sec. IIII B[ the ingoing fluxes 
do not affect the final spin, which only depends on the energy 
and angular momentum that leave the binary system. How- 
ever, the energy and angular momentum fluxes falling into the 
horizon might in principle spin the black hole up and even 
create a naked singularity before the particle is captured. 
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Because our code only calculates the outgoing fluxes, we 
cannot test this conjecture explicitly. However, we can notice 
that the analytic derivation of our scaling for the fluxes [Eqs. 
(fT9]l-(l22]i1 applies both to the outgoing and ingoing fluxes, 
because nowhere in the derivation we do make use of the 
fluxes being ingoing or outgoing. (Also, as already stressed, 
Ref. ll45ll found that the scaling ( fT6l ) from which we start holds 
not only for the outgoing fluxes but also for the total ones.) 
Because we have validated our scaling [Eqs. (fT9]l-(l22]i1 by 
comparing it to numerical results for the outgoing fluxes, we 
can now use it to assess the effect of the ingoing ones. In par- 
ticular, the spin of the black hole before the particle is captured 
is given by 

1 - 2e2 + 2V3e£;i"Afcycics , (36) 

where we have used Eqs. (I20li- (l22] | to estimate the ingoing 
fluxes E'lad.in and Liad.in- Using expressions ( fT9] l for and 
([T5ll for A^cycios, we obtain Li-ad.in - 2£'rad,in ^ log(fce). 
As for the outgoing fluxes, this scaling shows that for fixed 
k the ingoing fluxes are negligible with respect to the term 
2e^ in Eq. ( |36] |. if the initial spin is sufficiently large. In 
other words, if the initial spin is larger than some critical value 
a^^ij, no naked singularity can form before the particle is cap- 
tured. However, if a < a^^j^, the effect of the ingoing fluxes 
is dominant over the term quadratic in e. The exact value of 
a^j.ij will depend on the normalization of the ingoing fluxes 
^-rad.in and irad.in- Howcver, sincc the ingoing fluxes are ex- 
pected to be comparable to the outgoing ones - because the 
fluxes for JS orbits are produced when the particle sits at the 
light-ring, which roughly corresponds to the maximum of the 
effective potential for gravitational waves - a^j.;j should be 
similar to the critical spin acrit ~ 0.9998 (corresponding to 
e w 0.01 H^) relevant for the outgoing fluxes. 

Even more worrisome, the ingoing energy flux must 
be positive. This is because the ingoing energy flux can be 
negative only in the superradiant regime, which would require 
that the particle's orbital frequency be smaller than the hori- 
zon frequency. Since 17ph w 1/2 — (■\/3/2)e is larger than 
the horizon's frequency r^hor ~ 1/2 — e, it must be > 0, 
which implies that the intermediate spin (l36T l is larger than the 
initial spin. In particular, for a < a^^j^, the ingoing fluxes 
would overspin the black hole past the extremal limit and cre- 
ate a naked singularity before the particle is captured. We will 
see in the next section how the conservative self-force pro- 
vides a mechanism which has the correct order of magnitude 
and (possibly) the right sign to prevent both the formation of 
naked singularities due to the ingoing fluxes (for a < a'^^^^.), 
and the formation of naked singularities due to the particle's 
capture (for a > a'^^-^J. 

A. Suppressing dissipative effects with a ring of particles 

The above analysis shows that dissipative effects are not 
sufficiently strong to invalidate JS analysis. In fact, there is 
another simple argument indicating that dissipative effects can 



almost be neglected. Ours and JS's analysis considered a test 
particle of energy E, rest-mass fi and angular momentum L. 
Take now instead N particles each of energy E/N, rest-mass 
fi/N and angular momentum L/N. Throw these N particles 
separated by 2Tr/N radians along the equator JS calculation 
proceeds in exactly the same way, and the same results will 
be obtained. However, radiation effects will be suppressed. 
In fact, when iV — ^ oo, radiation should be suppressed com- 
pletely. The reason is that the ring can be thought of as a 
sum of point particles, and the fluxes result from the interfer- 
ence between the gravitational waveforms from all the parti- 
cles. Generically, this interference always reduces the energy 
output relatively to a single particle llsSmlll . 

In fact, for rings with angular momentum close to the criti- 
cal one, most of the radiation is emitted in a quasi-circular or- 
bit, as we have argued before. Because the contribution from 
particles in circular orbits has terms of the form S{uj — mil) 
and because m = for a ring, then the radiation in the circular 
regime is actually totally suppressed. This has been verified 
numerically by a number of authors. For trajectories plunging 
into a Schwarzschild black hole Oohara and Nakamura found 
that radiation is suppressed as the angular momentum of the 
ring is increased fell], while Kojima and Nakamura found sim- 
ilar results for plunging rings into rotating black holes ll60ll . 

Thus, the construction of a ring of JS-like particles is able 
to suppress radiation, while violating the Cosmic Censorship 
Conjecture, just like in the original JS analysis. This is our fi- 
nal argument showing why dissipative effects cannot prevent 
naked singularities from forming in this scenario, and that in- 
stead conservative effects must be taken into account. 



V. THE CONSERVATIVE SELF-FORCE 

In the previous section, we have shown that the radiation 
reaction can prevent the formation of naked singularities only 
for some of the JS orbits. More specifically, the final spin a/, 
which includes the effect of the radiated energy and angular 
momentum, is smaller than 1 only for orbits with impact pa- 
rameter b extremely close to the light-ring's impact parameter 
6ph, and the difference \b — 6ph| needed to ensure a/ < 1 be- 
comes smaller and smaller as the initial spin a approaches 1, 
i.e. the radiation reaction becomes less and less effective near 
the extremal limit. 

Also, we have shown that in the cases in which a/ < 1, 
the gravitational-wave fluxes that enter the horizon before the 
particle is captured tend to spin the black hole up past the ex- 
tremal limit, forming a naked singularity. While the capture of 
the particle would then offset this spin-up and produce a final 
spin Gf < 1, it is unclear whether the capture would happen 
at all, because the black-hole horizon would have disappeared 
and the particle may simply be scattered by the naked singu- 
larity. Also, even if the particle were captured, a naked singu- 
larity would exist for a finite amount of time, and the Cosmic 
Censorship Conjecture would still be violated. 

We will now use the results of the previous section to esti- 
mate another effect, the so-called conservative self-force lf52l 
I62I - I65I1 . and show that its order of magnitude is sufficient 
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to solve both these problems and prevent the formation of 
naked singularities for all of the JS orbits, provided that it 
carries a certain sign. Because our approach can only esti- 
mate the order of magnitude of the conservative self-force, 
but not its sign, our result can be viewed as putting forward a 
simple test of the Cosmic Censorship Conjecture using self- 
force codes. Unfortunately, to date none of the existing self- 
force codes ll66ll can handle orbits that, like the JS orbits, 
are relativistic and around almost-extremal Kerr black holes. 
As we will show, however, the existing calculations for non- 
relativistic orbits around non-spinning black holes seem to 
hint at a conservative self-force sign consistent with no for- 
mation of naked singularities. While more general self-force 
codes than those available today will be needed to say the last 
word on whether the JS orbits might create naked singular- 
ities, this result seems to hint at the conservative self-force 
playing a crucial role in enforcing the Cosmic Censorship 
Conjecture. 

In order to introduce the concept of conservative self- 
force, let us consider a black hole with gravitational radius 
Rg = 2G^/c^, moving in a curved background spacetime 
with C ^ Rg? In order to study the motion of this black 
hole in a completely rigorous way, one would need to set up 
a proper initial value formulation, but a reasonable alternative 
for practical purposes is to use a matched asymptotic expan- 
sion ll62l - l64ll . In particular, near the black hole {i.e. at dis- 
tances r to the black hole smaller than some limiting radius 
r,; ^ £), the metric can be written as 

^internal = ,9bh + Hi{r / C) + H^ir j + ... , (37) 

where pbh is the metric of an isolated black hole and 
Hi{r/C), H2{r/CY are corrections due to the presence of 
the "external" background. 

Far from the black hole {i.e. for r > re, being a suitable 
radius ^ Rg), the geometry is that of the background space- 
time plus perturbations due to the black hole's presence, and 
the metric can therefore be written as 

^external — ^/background ^1 {Rg/C)+h2{Rg/Cf + ..., (38) 

hi and h2 being functions (of time and position) representing 
the perturbations produced by the black hole. 

Because Rg ^ C, there exists a region < r < ri where 
both pictures are valid and the two metrics can be matched. 
Doing so, one obtains that the black-hole equations of motion 
are ll62l - i65ll 

u^V^u^ = + /d^i,, + 0{RglCf , (39) 

where V is the Levi-Civita connection of the background 
spacetime. The terms f^^^^ ^^'^ /diss ^i^g/^)^ ^re 
dubbed the conservative and dissipative self-force. Remark- 
ably, Eq. (3% turns out to be the geodesic equation of a parti- 
cle in a "perturbed" metric g = g + h^, where is a smooth 



^ Tliis discussion is completely general because the motion of a black hole 
is the same as that of a particle with the same mass ^, at leading and next- 
to-leading order in Rg/ C |63| . 



tensor field of order 0{Rg/C): 

u^V^.u'' = . (40) 

In this equation, the Levi-Civita connection V and the 4- 
velocity are defined with respect to the "perturbed" metric 
.9 = 5 + f^^^ ™d can be interpreted as the (regularized) 
metric perturbation produced by the black hole itself. 

The dissipative self-force can be shown ll67ll to be equiva- 
lent to the effect of the energy and angular-momentum fluxes 
on the particle's trajectory, which we considered in the previ- 
ous sections. We can therefore use the results of Sec.llVlto say 
something about the scaling of the self-force. Considering for 
instance the energy lost in gravitational waves, from the defi- 
nition of the particle's energy E = ~pt and from Eq. ( [39] l one 
obtains 

dE/dr = -/i/f = 0{Rg/Cf . (41) 

Assuming now that the background spacetime is a black hole 
with mass M ^ £ ^ Rg, and specializing to orbits near the 
horizon, from the geodesies equation one immediately obtains 
dt/dr ^ rn /{r — rji) ll4lll . which implies 

dE/dt - (r - rH)0{RglCf . (42) 

Because in the extremal Kerr geometry the ISCO, the 
marginally bound circular orbit and the light ring coincide 
with the horizon [46], this equation will be valid for the circu- 
lar orbits between the ISCO and the light-ring in the almost- 
extremal case. Also, because for these orbits the energy E 
is proportional to the mass ji, we can think of Rg as being 
proportional to the energy E in Eq. ( |42] | (this corresponds to 
replacing /i with E when going from Eq. ( fTSI l to Eq. ([TtI i in 
Sec. 1111 Al) . We can then follow the same procedure that we 
used in Sec. 1111 Al i.e. we can consider a sequence of unsta- 
ble circular orbits between the marginally bound orbit and the 
light-ring, with constant positive energy (and therefore with 
rest-mass going to zero as the light-ring is approached). Tak- 
ing the limit of Eq. (|42] | along this sequence of orbits shows 
that Eq. ( |42] | should also be valid for a photon at the light-ring, 
provided that Rg is interpreted as proportional to the photon 
energy. 

As already stressed, the numerical results of Sec. |IV] sup- 
port this interpretation of Eq. (l42t . In fact, for the JS or- 
bits Rg should scale like the particle energy E ^ e, and 
because — ^ C'(e) we have that Eq. ( l42b gives 
dE/dt ~ 0{e)^, which is indeed the scaling that we found 
numerically in Sec.|lV] It therefore seems that for a black hole 
with E ^ fi, the size entering the matched asymptotic expan- 
sion analysis that we sketched above is Rg ~ 2GE/c? and 
not Rg = IG^ijc?. This is hardly surprising, since the size 
associated with a black hole or particle moving at relativis- 
tic speeds is given by its energy and not by its mass, simply 
because in General Relativity energy gravitates. 

Further evidence comes from the so-called Aichelburg-Sexl 
metric, which represents a Schwarzschild black hole as seen 
by an observer moving at nearly the speed of light. (More pre- 
cisely, the Aichelburg-Sexl metric can be obtained by boost- 
ing the Schwarzschild metric to the speed of light, keeping 
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the total energy fixed.) As one would expect from physical 
intuition, this metric depends on the total energy E and not on 
the rest-mass ll68ll . and in particular this boosted black hole 
absorbs particles within a distance ^ E from it. 

To make the argument more rigorous, one may attempt 
to set up a matched asymptotic expansion for a photon in a 
generic curved background. Because the stress-energy tensor 
of a photon [Eq. ( fTSl lI depends on its energy and not on its rest- 
mass (which is of course zero), and because the metric near a 
photon is presumably given by the Aichelburg-Sexl metric, 
the only size entering the matched asymptotic expansion and 
therefore the self-force will be i?g = 2GE/c^. Developing 
such a formalism goes beyond the scope of this paper, but the 
picture that we described above is physically clear, and it is re- 
markable that we were able to test it with the numerical results 
presented in Sec. |IV] 

Based on this argument, the size of a black hole or particle 
moving on a JS orbit is given hy E ^ e, which is sufficient to 
conclude that the conservative self-force affects JS's analysis. 
This is easy to see from Eq. (|40] | (althou gh t he same result 
can be obtained from Eq. ( [39] i: see Ref. ll69ll ): because the 
regularized metric "perturbation" is 0{Rg/£) ~ 0{e), 
the effective potential for the radial motion differs from the 
"geodetic" one by 0{Rg/C) = 0(e) iHsIItI. In particu- 
lar, the light ring's impact parameter &ph becomes &ph + Sbph, 
with Sbpii = 0{Rg/ C) = 0{e). Because the JS orbits have 
feph — b = 0{Rg/C) = 0{e), the conservative self-force may 
prevent them from plunging into the horizon. This effect is 
intuitive: if the particle's size is ^ e, finite-size effects are 
important for impact parameters b — bph + 0{e) (i.e. the im- 
pact parameter is so close to the light ring's impact parameter, 
which discriminates between plunging and scattering orbits, 
that finite-size effects must be taken into account). 

Clearly, what remains to be determined in this analysis is 
the coefficient of the impact parameter change (55ph produced 
by the conservative self-force. In particular, the sign of (5feph is 
crucial, because if dbph > the light ring's impact parameter 
would increase, actually making it easier for the JS orbits to 
plunge into the black hole and produce a naked singularity. 
If instead (56ph < 0, the light ring's impact parameter and 
therefore the black-hole photon cross section would shrink, 
thus making it harder for the JS orbits to hit the black hole. 

As we have already mentioned, a calculation of (5feph is 
not feasible with present self-force codes ||66l]. which can 
only handle non-relativistic orbits around non-spinning black 
holes. Nevertheless, we can try to use the existing results for 
these orbits to guess how the impact parameter of the light 
ring might change under the effect of the conservative self- 
force. In particular, Ref. j69| showed that the ISCO frequency 
in a Schwarzschild black hole increases due to the conserva- 
tive self-force. If the same behavior applies to relativistic or- 
bits in almost extremal Kerr black holes, i.e. if fiph for these 
spacetimes increases under the effect of the conservative self- 
force, then 6pii should decrease {i.e., (56ph < 0), because for 
circular photon orbits one has 6ph = 1/flph. (This follows 
from the fact that the photon 4-momentum is a null vector: 
from p'^''p^ = and from p^ = —E5j^ + LSf^ and p'' — 
p^{S^ + fiiJ^), one immediately obtains b ^ L/E = 1/ft.) 



This would imply, as already mentioned, that the black-hole 
photon cross section would shrink under the effect of the con- 
servative self-force, possibly preventing the JS orbits from be- 
ing captured and naked singularities from being formed. 

Clearly, explicit self-force calculations for relativistic or- 
bits around almost extremal Kerr black holes will be needed 
to confirm this conjecture and determine the exact numerical 
value of (56pii (in order to go beyond the order of magnitude re- 
sult Sbph = 0(e) derived above). It is of course very well pos- 
sible that such explicit calculations will find that 6ph increases 
rather than decreases for almost extremal black holes, or that 
its decrease is too small to prevent the JS orbits from being 
captured. Hints at a possible change of sign of Sbph when go- 
ing from a = to a ~ 1 come for instance from Ref. iItTI . in 
which Warburton and Barack found that the ISCO frequency 
for a scalar particle increases under the effect of the (scalar) 
self-force when \a\ < 0.9, but decreases for a > 0.9. A sim- 
ilar result was found by Refs. ll72ll . who calculated the ISCO 
shift due to the conservative self-force usin g an effective one- 
body model for spinning black-holes 11721 \TM . calibrated in 
the a = case with the results of Ref. 116911 . and found that the 
ISCO frequency decreases at high spins. These results, albeit 
still inconclusive, highlight even more compellingly the need 
for a rigorous calculation of the gravitational self-force for ul- 
trarelativistic particles (or photons) moving in an almost ex- 
tremal Kerr background. Such a calculation will probably be 
a very significant step towards an understanding of the range 
of validity of the Cosmic Censorship Conjecture. 

If, however, we assume that (56ph decreases in the almost- 
extremal limit under the effect of the gravitational self-force, 
a change of order (56ph = 0{t) may also be enough to prevent 
the ingoing fluxes from falling into the horizon and create a 
naked singularity. Suppose in fact that 8bp\^ were negative 
enough to prevent ultrarelativistic particles with impact pa- 
rameter b in the JS range (6 = 2 + (5e with 2\/2 < 5 < 4, 
cf. Sec. from falling into the horizon. These particles will 
then describe a large number of quasi-circular orbits near the 
periastron before being scattered. During these quasi-circular 
orbits, they will emit gravitational fluxes E'lad and Ljad into 
the horizon, with L^^^jE^^^ = l/ilpdi = b, and these fluxes 
tend to spin the black hole up, possibly above the Kerr limit 
\cf. discussion around Eq. ([36ll1. However, it is well-known 
that in the eikonal limit {i.e., in the small-wavelength limit, 
corresponding to large m's in the decomposition of Sees. 1111 B] 
and 1111 Cl l gravitational waves behave like massless particles 
("gravitons"), their propagation in a Kerr spacetime being 
regulated by a radial effective potential that is the same as 
the radial effective potential regulating the motion of pho- 
tons [c/ Eq. (O].** We can therefore think of the ingoing 
fluxes as being made of wavepackets with impact parameter 
fog ranging from to +oo, but with a distribution centered 
(roughly) on the impact parameter b of the particle produc- 



* The eikonal approximation is a very good one for the gravitational waves 
emitted by the JS orbits because, as we stressed in Sec. IIII CI the fluxes 
in the large-m multipole moments turn out to be very important for these 
orbits. 
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ing them (because the ratio of the total ingoing fluxes must 
be Liad/^'rad = b). Clearly, the only wavepackets capable of 
overspinning the black hole are those with bg in the JS range 
bg = 2 + Se with 2a/2 < (5 < 4. However, if the effective po- 
tential for ultrarelativistic orbits (and therefore for gravitons) 
changes under the effect of the self-force in such a way as to 
prevent the JS orbits from being captured [i.e. if Sbp^ ~ Oit) 
is sufficiently negative], it will also scatter the gravitational 
wavepackets with hg in the JS range. 

We stress that Jacobson and Sotiriou proposed also a dif- 
ferent scenario in which naked singularities might form. In 
particular, they consider a spinning particle having energy E 
and spin S, with the energy satisfying the limits (|2]i and the 
spin satisfying the limits ([T]i (where we identify L with the 
particle's spin S). They also assume that the particle has spin 
parallel to the black-hole spin and that it is dropped into the 
black-hole horizon along the common direction of their spins. 
The final spin of the resulting black hole will be larger 
than 1, because the lower bound in Eq. ([T]i is achieved pre- 
cisely by imposing a?^ > 1. Also, the condition for a spin- 
ning particle to fall into a Kerr black hole is i? > flnS IBol . 
where Hh ~ 1/2 — e is the horizon frequency, and the up- 
per bound of Eq. ([T]i can indeed be written in that form, i.e. 
the upper bound of Eq. ([U ensures that the spinning particle 
actually falls into the black-hole horizon. 

From the discussion above, however, it is clear that the 
conservative self-force will modify the metric and therefore 
the horizon frequency by terms of order E ^ e. In par- 
ticular, if the horizon frequency increases under the effect 
of the conservative self-force (just like the ISCO frequency 
increases in a Schwarzschild spacetime lHH) and becomes 
Hh = + K£ ~ 1/2 — (1 — K)e (with k > being a 
coefficient), Eq. ([T]i would become 

2e^ + 2E + E^ <S <[2 + A{l-K)e]E. (43) 

Imposing that the upper bound in this equation be larger than 
the lower bound [i.e., imposing that the interval described by 
Eq. ( |43] l is not empty], we immediately obtain that if k > 
1 - I/V2 w 0.293 there are no orbits that both fall into the 
black hole and create a naked singularity. In other words, the 
conservative self-force may be enough to prevent the JS from 
forming naked singularities even in the case of particles with 
spin. 

Finally, proceeding exactly in the same way as we just did 
for spinning particles, it possible to show that the self-force 
may also be enough to prevent naked singularities from form- 
ing in the case of non-spinning particles on bound JS orbits. 



As we stressed in Sec. [Ill these orbits start from very close to 
the black-hole horizon and do not orbit the light ring multiple 
times, so the analysis of the previous sections does not apply 
to them. However, because their energy and angular momen- 
tum must still satisfy Eqs. ([T]i and (|2|i, the horizon frequency 
will change by 0{e) due to the conservative self-force. In 
particular, if the horizon's frequency increases, then the up- 
per bound in Eq. ^ will be lowered and the allowed angular 
momentum interval might possibly shrink to nothing. 
VI. CONCLUSIONS 

As discussed in the introduction, there is strong circumstan- 
tial evidence for the stability of Kerr black holes, which leads 
us to believe that rumors of their death may have again been 
greatly exaggerated. Nevertheless, the particular mechanism 
proposed by Jacobson and Sotiriou ||36|] is exciting enough to 
deserve careful thought. The mechanism can be at play in as- 
trophysical settings and the understanding of why it fails (or 
not) will certainly shed light on highly dynamical processes 
close to extremal Kerr black holes. We have shown that radia- 
tion by point particles close to the last stable circular geodesic 
is suppressed as the black hole approaches extremality, pro- 
viding further support to earlier results by Chrzanowski and 
Misner ll42ll43ll . Indeed, we have shown that if one modifies 
Jacobson and Sotiriou's analysis by replacing the test parti- 
cle with a ring of particles, the gravitational radiation is sup- 
pressed by interference effects. It therefore seems that the 
conservative self force might be the main effect preventing vi- 
olations of the Cosmic Censorship Conjecture. While we have 
provided arguments in favor of this picture, a rigorous proof 
is still unavailable, and the role of the conservative self-force 
in the Cosmic Censorship Conjecture remains an outstanding 
open issue. 
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